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Abstract 

Observation of the neutrinoless double beta decay (Ov/3/3) has es- 
tablished that there is lepton number violation in nature and the neu- 
trino masses are Majorana in nature. It also gives the absolute mass 
of the neutrinos and discriminates between different models of neu- 
trino masses. The allowed amount of lepton number violation puts 
severe constraints on some possible new physics beyond the standard 
model. The recent results from WMAP are consistent with the con- 
sequences of the neutrinoless double beta decay. They improve some 
of these constraints very marginally, which we shall summarise here. 
We mention the new physics which are not affected by WMAP and 
could make both these limits from the neutrinoless double beta decay 
and WMAP consistent. 



1 Introduction 



During the past few years there have been several new results in neutrino 
physics [1-6]. The atmospheric neutrino problem established that there is 
neutrino mass and that oscillations occur between — > v T . This is sup- 
ported by K2K [1]. The solar neutrino problem started quite some time 
back and experiments favored the large mixing angle MSW solution [2] . The 
KamLAND experiment has now confirmed that the large mixing angle MSW 
solution is the solution of the solar neutrino problem [3] . These experiments 
determine the two mass-squared differences and two mixing angles of the 
neutrino mass matrix. There is only an upper bound on the third mixing 
angle coming from the reactor experiments [4] . The neutrinoless double beta 
decay tells us that the neutrinos are Majorana particles and also provide us 
with the absolute mass [5]. Recently the Wilkinson Microwave Anisotropy 
Probe (WMAP) has provided us with a value for the total mass of the neu- 
trinos and claims that there are three generations of neutrinos [6]. We shall 
restrict ourselves to only three generations of neutrinos. 

A fit of all the data for atmospheric neutrinos and the K2K gives [7] 

Ami - (1-8 - 4.0) x 10" 3 eV 2 , sin 2 29 atm > 0.87. (1) 

This oscillation is established to be between a and a v T and the possibility 
of a sterile neutrino is ruled out. A global fit to the results from the atmo- 
spheric, solar and the reactor neutrinos gives us two solutions for the solar 
neutrinos [8, 9]. The allowed 3a region for the two solutions are 

LMA-I solution: 5.1 x 10~ 5 eV 2 < Am 2 oil < 9.7 x 10~ 5 eV 2 

LMA-II solution: 1.2 x 10~ 4 eV 2 < Am 2 sol2 < 1.9 x 10" 4 eV 2 
with mixing angle 

0.29 < tan 2 9 sol < 0.86. 

For increasing sin 2 #13 the LMA solutions shrink and eventually disappear, 
which requires sin 2 9 13 < 0.04 [8]. Again, a sterile neutrino solution is com- 
pletely ruled for the solar neutrinos. 

An analysis of the Heidelberg-Moscow data yields a half-life for the neu- 
trinoless double beta decay experiment of [5] 

T°Y 2 = (0.7 - 18.3) x 10 25 y at 95 % c.l. 
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with a best value of 1.5 x 10 25 y. The signal for the neutrinoless double 
beta decay amounts to an effective Majorana neutrino mass of the electron 
neutrino in the range of [5] 

(m) = (0.11 - 0.56) eV at 95% c.l. (2) 

with a best value of 0.39 eV (with nuclear matrix element of reference [10]). 
A weaker bound for the neutrino mass of 

(m) = (0.05 - 0.86) eV at 95% c.l. (3) 

is deduced, when a ±50% uncertainty in this nuclear matrix element is al- 
lowed (for details see [5]). We shall not discuss the bounds on the neutrinoless 
double beta decay from another analysis [11] that uses a 15-years old nuclear 
matrix element. In that calculation of the nuclear matrix elements they did 
not include a realistic nucleon-nucleon interaction, which has been included 
by all other calculations of the nuclear matrix elements over the last 15 years. 
This matrix element has been ruled out by the WMAP result completely [12]. 

On the basis of the most recent result from WMAP it is claimed that 
there are three generations of neutrinos and that the LSND result is ruled 
out [12]. A limit on the total neutrino masses of 

m s = J2m u < 0.69 eV at 95% c.l., (4) 

is given by the analysis of ref. [6]. It has been shown, however, that this 
limit may not be very realistic. Another analysis shows that this limit on the 
total mass should be [13] 

m s = Y j m v < 1.0 eV at 95% c.l. (5) 

The latter analysis also shows, however, that four generations of neutrinos 
are still allowed and in the case of four generations the limit on the total mass 
is increased to 1.38 eV. If there is a fourth neutrino with very small mass, 
then the limit on the total mass of the three neutrinos is further weakened 
and there is essentially no constraint on the neutrino masses. In our analysis 
we comment on these possible values. 
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2 Constraints on models of neutrino masses 



There are several consequences of the neutrinoless double beta decay. Fol- 
lowing the announcement of the positive evidence for the process, lots of 
activity started (see e.g. [14]). We concentrate on a few aspects only. We 
first discuss the different scenarios of neutrino masses. We parametrize the 
neutrino mass eigenvalues as 

A ™L = l A ™i2l and Am 2 atm = \Am 2 23 \ (6) 

where Am\ 2 = \rri2\ 2 — \mi\ 2 and Am 23 = m\ — \ rri2\ 2 . For the mixing angle in 
the atmospheric neutrinos we assume maximal mixing, so that sin 2 29 a t m = 
sin 2 2^23 = 1, sin ^23 = cos #23 = l/\/2- F° r solar neutrinos we allow the 
entire allowed range and write cos# 12 = c and sin# 12 = s, with 0.54 < s/c < 
0.93 or 0.48 < s < 0.68 and a best value of s = 0.55. Taking the limit on 
the third mixing angle, we assume, cos #13 ~ 1 and sin #13 — u < 0.2. In our 
analysis we do not include CP violating phases and hence the mixing matrix 
may be parametrized as 



c s u ^ 

U= I -(s + cu)/V2 (c-su)/V2 l/y/2 
(s-cu)/V2 -(c + su)/V2 l/y/2 j 



(7) 



As noted in earlier references [15], inclusion of CP violation does not affect 
most of the conclusions regarding the bounds on the neutrino mass. The 
neutrinoless double beta decay bound then implies 

mic 2 + m 2 s 2 + m 3 u 2 = (m), (8) 

and the WMAP implies 

|^i| + \m 2 \ + \m 3 \ < m s . (9) 

In some models it is possible to have neutrinoless double beta decay mediated 
not be exchange of a massive Majorana neutrino, but of some exotic particles 
[16, 17]. In that case it would make no sense to compare limits from WMAP 
or corresponding experiment with neutrinoless double beta decay. 

In general, it may be possible to classify the neutrino masses in four 
classes: 
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Hierarchical: 

This is the most natural choice for neutrino masses, where all three neu- 
trino masses are different and hierarchical (m 1 < m 2 < to 3 ) similar to the 

mass hierarchy of the charged fermions. This implies m 3 = m atm = \J Aml tm , 

m 2 = m so i = \J Am 2 sol and m 1 < m so [. The WMAP bound then becomes, 
1^3 1 < fn S i since the other two masses are too small. If we use the highest 
value of mi to be m so i and use the WMAP constraint along with CHOOZ and 
solar neutrino results, then the contribution to the effective mass appearing 
in the neutrinoless double beta decay becomes 

(to) < m so i + m 3 u 2 < m so i + m s u 2 < 0.027 eV 

for m s = 0.69. Thus the WMAP constraint sets a stronger limit compared to 
the neutrinoless double beta decay. But it still allows the hierarchical solu- 
tion with solar and atmospheric neutrino solutions, whereas the neutrinoless 
double beta decay result rules out this solution. This is true also for a higher 
value of the WMAP constraint of m s = 1.0, where (to) < 0.04 eV. How- 
ever, when the atmospheric neutrino constraint is considered m 3 = m atm , the 
WMAP condition is trivially satisfied. In this case the contribution to the 
neutrinoless double beta decay is further reduced and becomes too small to 
explain the present bounds (to) < m so i + m atm u 2 < 0.016 eV. We considered 
the maximum value for mi ~ m so \ and the largest values of all the parame- 
ters. There is no lower limit on the contribution to (to) in this case. In brief, 
neutrinoless double beta decay does not allow the hierarchical neutrino mass 
matrix. The WMAP limit cannot yield such a strong statement. 
Degenerate: 

This is the most interesting solution at present, which is allowed by all 
the experiments. Here one assumes m\ m 2 ~ m 3 w m , where m is the 
overall mass and the mass squared differences are as required by the solar and 
atmospheric neutrinos. The WMAP bound then implies m < to s /3. The 
neutrinoless double beta decay constraint now implies (to) < m < m s /3 
and hence part of the allowed region is ruled out by the WMAP result. For 
m s = 0.69 eV, this implies (to) < 0.23 eV and the best fit value for the 
neutrinoless double beta decay is not within the range. On the other hand, 
for the other more realistic limits of m s = 1.0 eV and 1.38 eV, we get (to) < 
0.33 eV and 0.46 eV respectively, which does not conflict with the best fit 
value of the neutrinoless double beta decay. For the degenerate solution the 
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lowest contribution to the neutrinoless double beta decay would correspond 
to m = m atm and comes out to be (to) > m atm (c 2 — s 2 — u 2 ) > 0.001 eV. 
Inverted Hierarchical: 

In this case one considers that two of the neutrinos are degenerate and 
heavier than the third one, toi ~ to 2 > to 3 . The Ato 2 2 is very small and 
can explain the solar neutrino problem, while Toi ~ vn atm so that Ato 23 can 
explain the atmospheric neutrino problem. The WMAP bound now implies 
that Toi = (toi + to 2 )/2 < m s /2 = 0.345 eV for m s = 0.69. In this case the 
effective mass for the neutrinoless double beta decay becomes 

(to) = toic 2 + to 2 s 2 < toi < m s /2 < 0.345 eV. 

for to s = 0.69. For m s = 1.0 the bound is 0.5 eV. The WMAP result does 
not allow part of the allowed region of the neutrinoless double beta decay, 
which is anyway not allowed by the atmospheric and solar neutrinos. The 
atmospheric neutrino solution requires Toi ~ m atm , so the contribution to 
the neutrinoless double beta decay becomes (to) < toi ~ m atrn < 0.06 eV. 
This value is still marginally allowed, when allowing for a ±50% uncertainty 
in the matrix element of [5]. The lowest allowed value now corresponds to 
(to) > m x {c 2 - s 2 ) = 0.003 eV. 
Partially Degenerate: 

This scenario is ruled out by neutrinoless double beta decay, but still 
allowed by atmospheric and solar neutrinos and also WMAP. There are 
two degenerate neutrinos whose mass difference squared Ato 2 2 solve the so- 
lar neutrino problem and the third neutrino is heavier to x w to 2 < to 3 . 
A solution to the atmospheric neutrinos requires to 3 = m atm . The main 
difference between the hierarchical and this scenario is that, in this case 
toi rs to 2 > m so \. To distinguish this solution from the degenerate sce- 
nario, we consider toi < m atm = 0.4 eV. The WMAP constraint now implies 
to 3 < to s = 0.69 eV. In this case we can have Toi ~ to 2 to be close to but 
less than to 3 < m s . So, the effective mass for the neutrinoless double beta 
decay becomes 

(to) = toic 2 + m 2 s 2 + m s u 2 < toi < m s = 0.69 eV. 

For m s = 1.0 this bound is also high as shown in the figure. If we now include 
the atmospheric neutrino result, the WMAP constraint is satisfied. The neu- 
trinoless double beta decay contribution now becomes (to) = toic 2 + to 2 s 2 < 
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mi < m atm - The highest value of TOi and TO2 correspond to to-i m m atm , 
which is included in the degenerate solution. So, if we restrict the partially 
degenerate solution to m\ < m atm , then the contribution to the neutrino less 
double beta decay becomes (to) < 0.04 eV and hence the neutrinoless double 
beta decay does not allow this solution, even when we relax the uncertainty 
in the nuclear matrix elements. 

All the four scenarios we mentioned are solutions of the solar and at- 
mospheric neutrino problems. All these solutions are also allowed by the 
WMAP result, except in the degenerate case when a small part of the solu- 
tion is not allowed. On the other hand, the neutrinoless double beta decay 
does not allow the hierarchical and the partially degenerate solutions. The 
inverted hierarchical solution is only marginally allowed by the neutrinoless 
double beta decay when an extra ±50% uncertainty in the nuclear matrix 
element is permitted. Only the degenerate solution is allowed by the entire 
neutrinoless double beta decay range, but a small part of the allowed region 
is ruled out by the WMAP result. This is shown in figure 1. In the figure we 
considered the more recent analysis of the WMAP result [13] and used the 
bound Tfi s < 1.0 eV. 

3 Constraints on lepton number violation 

The neutrinoless double beta decay process could be triggered also by ex- 
change of other particles, than massive neutrinos [21]. In this sense a de- 
duced effective mass is - though the most natural explanation - strictly only 
an upper limit. The measured half-life (or its lower bound) can thus be used 
to deduce limits for other beyond standard model physics and other lepton 
number violating interactions. 

The WMAP constraint depends on which value we consider. If we con- 
sider the value quoted in the original paper, then it implies a slightly im- 
proved bound on some of the lepton number violating processes. However, 
if the weaker bounds are considered, then it does not improve any of the 
constraints compared to the present bounds on the neutrinoless double beta 
decay. In the rest of the analysis we shall consider the value J2 m iy < 1-0 eV, 
so that combining with the neutrinoless double beta decay we get a limit on 
the effective mass of (to) < 0.33 eV. 

If there are heavy right-handed neutrinos, which have small mixing with 
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Figure 1: Contributions in different models to the neutrinoless double beta 
decay. The present result is given by the dark shaded region (the solid line de- 
noting the best value and the light shaded region allowing ±50% uncertainty 
in the nuclear matrix element). The WMAP line is plotted for J2 m v < 1-0 
eV, although it could be even weaker (i.e., the line running higher) as men- 
tioned in the text. Future sensitivity that might be reached for the CUORE 
[18], MOON [19] and the one ton and ten tons GENIUS [20] are given for 
comparison. 



the left-handed neutrinos, then they can enter the neutrinoless double beta 
decay processes and contribute to the effective mass. The present bound on 
the lifetime of the neutrinoless double beta decay would then give a constraint 
[22] 

M N >6x 10 7 GeV. (10) 

From the same analysis the bound on the the right-handed W boson comes 
out to be 
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With some reasonable theoretical input, this may be translated to an absolute 
lower bound of mw R > 1-2 TeV [22]. Including WMAP constraints this 
bound will be improved to iti Wr > 1.5 TeV. 

Using the lifetime of neutrinoless double beta decay, the probability for 
the discovery of the inverse beta decay process e~e~ — > W~W~ at NLC 
could be constrained. The present value [5] can be achieved at a future linear 
collider NLC when it reaches a center of mass energy of 2 TeV [23, 24]. There 
is hardly any change in the analysis when the WMAP result is considered. 

If there are Higgs scalar bilinears, which couple to the usual quarks and 
leptons, they can also allow for neutrinoless double beta decay. However, in 
this case it is possible that the contribution of these scalars to the neutrino 
mass is negligible [16]. Then these scalars could allow the neutrinoless double 
beta decay as claimed, but the neutrino mass will be much smaller and hence 
there will not be any constraint from the WMAP result. In other words, in 
these scenarios it is possible to satisfy even the stronger WMAP bounds 
simultaneously making it consistent with the neutrinoless double beta decay. 

Although the dileptons were first considered [25] in connection with the 
left-right symmetric model, no significant bound is possible on this scalar. 
For the leptoquarks the bound is better. If the X— type leptoquarks (SU (2) l 
singlets) mix with the Y — type leptoquarks (SU(2)l doublets), then they can 
give an effective operator uvdl that generates a diagram contributing to the 
neutrinoless double beta decay involving the leptoquarks [26]. A mixing 
between these two leptoquarks could take place only after the electroweak 
symmetry breaking, if both these leptoquarks couple to the usual standard 
model Higgs doublet <p. In that coupling (j>XY will induce a mixing of 

X with Y when <fi acquires a vev. It was noticed [26] that in the leptoquark 
mediated case, there is a huge enhancement factor of ~ 10 s (leV/m„), 
where m u is the effective neutrino mass entering the neutrinoless double beta 
decay contribution, and < q > is the Fermi momentum of a nucleon inside 
a nucleus, which is about 200 — 300 MeV. For a leptoquark with mass of 
the order of 100 GeV, the effective coupling constant (including the mixing 
contribution) comes out to be about 10 -9 . The WMAP result does not 
modify this bound. 

There are other exotic scalar bilinears, which can also mediate the neu- 
trinoless double beta decay [27]. From the present allowed range of the 
neutrinoless double beta decay it is thus possible to put severe constraints 
on these scalars. If we assume a common mass for these scalars of about 
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100 GeV and that the self interaction of these scalars is of the order of 1, 
a strong bound on the effective coupling of the scalars to ordinary fermions 
becomes / < 10~ 7 . In other words, if the couplings are assumed to be of 
the order of 1, there is a lower bound on the masses of these scalars to be of 
the order of 10 4 GeV. In general, a constraint on the ratio of the masses to 
their coupling to first generation fermions for all the exotic scalar bilinears 
could be given [27] , which is comparable to the bounds from other processes 
like K° — K° oscillations, B° — B° mixing, D° — D° mixing, proton decay or 
n — n oscillations. Again, these bounds are not modified by the WMAP re- 
sult. Only in some cases, when the exotic particles also give a large neutrino 
mass, the present bound is improved very marginally by the WMAP result. 
Considering the uncertainty in the matrix elements for these processes, these 
bounds are not worth mentioning. 

The indirect bounds discussed also constrain some of the possibilities of 
composite particles. For example, if a neutrino is a composite particle, the 
most severe constraint comes from the neutrinoless double beta decay [28], 



where A c is the compositeness scale, M N is the mass of the heavy excited 
neutrino and / is the dimensionless coupling constant. This bound is also 
not affected by the WMAP result. 

It is also possible to constrain several parameters of supersymmetric the- 
ories. Recently one analysis claimed that the WMAP result improves the 
existing bounds on these parameters by about one order of magnitude [29], 
but we find this claim unreasonable. If we consider the bound from the 
neutrino mass from the neutrinoless double beta decay, then there is hardly 
any change in this bound after including the WMAP result. For the bounds 
available in the literature for the supersymmetric theories coming from the 
neutrinoless double beta decay see for example, ref. [30, 31]. Given the un- 
certainty in the calculations of these models, this change in the number is 
negligible. We do not present here these unchanged numbers. 

The bounds on the sneutrino-antisneutrino oscillation is also unchanged 
compared to the earlier bounds [32]. There are also bounds on the scale of 
extra dimensions in models in which mini black holes generate neutrino mass 
[33]. That limit is modified very marginally. The allowed textures of the 
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neutrino masses are also constrained by the neutrinoless double beta decay, 
which are not affected by the WMAP results [34]. 

4 Summary 

We studied the consequences of the neutrinoless double beta decay and the 
WMAP results. Models of neutrino masses are severely constrained by the 
neutrinoless double beta decay, while the WMAP result may only eliminate 
a small part of the region allowed by the neutrinoless double beta decay. 
WMAP constraints become the strongest in the case of hierarchical solu- 
tion, but it is not ruled out. On the other hand, the hierarchical solution is 
ruled out by the neutrinoless double beta decay. The degenerate solution is 
most favored by the neutrinoless double beta decay and only a small part 
of the allowed region is constrained by the WMAP result. WMAP bounds 
become weaker for the inverted hierarchical and partially degenerate solu- 
tions, whereas neutrinoless double beta decay can only marginally allow the 
inverted hierarchical solution when an extra ±50% uncertainty in the nuclear 
matrix element is allowed. The constraints on the lepton number violating 
processes are severely constrained by the neutrinoless double beta decay, 
while only a few of these constraints are negligibly modified by the WMAP 
result. 

Acknowledgement One of us (US) thanks the Max-Planck- Institut fur 
Kernphysik for hospitality. 
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